Gluon Condensate, Wilson Loops and Gauge/String Duality 
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We test gauge/string duality by evaluating expectation values of small Wilson loops in pure 
Yang-Mills theories. On the gauge-theory side, there exists rich phenomenology. The dual formula- 
tion provides with a universal language to evaluate the gluon condensate and quadratic correction 
in terms of the metric in the fifth coordinate. Quantitatively, the estimated value of the gluon 
condensate is approximately 0.010 GeV 4 . 
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PACS numbers: 11.15.Tk, 11.25Pm, 12.38Lg 

Introduction. The determination of the condensates 
of the QCD vacuum is a very important issue in 
phenomenology of strong interactions As known, 

the condensates can only be determined in a non- 
perturbative formulation of the theory. Until recently, 
the lattice formulation was in this respect a unique theo- 
retical tool. There is a long history of attempts to deter- 
mine the gluon condensate from first principles by com- 
puting the expectation value of the gluonic action and 
subtracting the perturbative contribution 0, S|. One 
needs, however, to subtract many orders of perturbation 
theory, and the interpretation becomes less straightfor- 
ward. Thus, there is a strong need for new alternative 
approaches to the problem. 

This subject has taken an interesting turn with Mal- 
dacena duality Q. Although the original proposal was 
for conformal theories, various modifications have been 
found that produce gauge/string duals with a mass gap, 
confinement, and supersymmetry breaking 

In this paper we will address the issues of the gluon 
condensates and of the quadratic correction in a dual 
formulation, known as AdS/QCD or holographic QCD. 

The model. Let us first explain the model to be con- 
sidered. The dual string spacetime is the product of five- 
dimensional space, with the Euclidean metric [6| 

ds 2 = c? cz (cbW + dz 2 ) , (1) 

and some five dimensional internal compact space X. We 
also take a constant dilaton. The AdS/QCD approach 
is a simplified version of duality that assumes a trivial 
dependence on the internal space X Q • Unlike other du- 
als, this model does share a few key features with QCD 
that singles it out and makes it very attractive for phe- 
nomenology: First, it is a nearly conformal theory at UV. 
Second, it results in linear Regge-like spectra for mesons 
[1, 0] . This fact allows one to fix the value of c. It is of 
order 0.9 GeV 2 Finally, the model results in a phe- 
nomenologically satisfactory description of the confining 
potential as well [ll| . Note that the model does not con- 
tain any free fit parameter. Thus, evaluation of the gluon 



condensate we are going to undertake can be considered 
as a next crucial test of the model. 

Given the metric, it is straightforward to calculate the 
expectation value of a Wilson loop by using the prescrip- 
tion of AdS/CFT 0. We first need to find the minimal 
area of string worldsheet which ends on a loop C living on 
the boundary (z — 0) of spacetime. Then the expectation 
value of the loop is given by 



(W(C))^e- s . (2) 
Finally, the gluon condensate introduced in [l[ 
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is given by the coefficient of s 2 , with s the area of the 
loop, in the expansion of the Wilson loop as s — > 0. Note 
that higher dimensional condensates are determined as 
coefficients of s n for n > 2. 

Calculating the circular Wilson loop. Now, let C be a 
circular loop of radius a. The worldsheet area is given by 
the Nambu-Goto action S — / d 2 £ ^y, with the 
induced metric on the worldsheet. We choose £ 1 = r and 
£ 2 = tp, where (r, ip) are the polar coordinates on a plane 
in R 4 . After introducing a new variable t = Wl — (^) 2 
and setting tp = (f ) 2 , the action becomes 



S = q I dttijj-2 e x ^\lil) + 



1-t 2 
At 2 



W) 2 , (4) 



where Q = and A = ^ca 2 . A prime denotes a deriva- 
tive with respect to t. The Euler-Lagrange equation com- 
ing from this action is 



1 / At 2 \ 

W> - - a^((v/) 2 + 7^2^) 



(•5) 



b) =0. 
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We are interested in the physical situation where the 
loop may be considered as relatively small, with a of or- 
der 0.05 fm. This means that the value of A is of order 
0.03. Thus, we can expand ip in powers A such that 
tp = J2 n =o VViA™, with the ip n 'B obeying a set of differen- 
tial equations. 

In particular, ipo is determined from 



This equation has the first integral 



= 0. 

(6) 
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We look for solutions of © which are regular and non- 
vanishing at t = 1. It follows from that 7 = on this 
class of solutions. Since tpo is positive for t < 1, the non- 
linear equation reduces to a simple linear differential 
equation 



2t - tp' = 



(8) 



whose solution with the boundary condition i/jq (0) = is 
simply 



(9) 



In the original vari- 
2 and describes the 



This solution was found in \12\, I X 3 1 
ables it takes the form z — \J a 2 — r 
minimal surface in AdSs bounded by a circle of radius a. 
Note that it can be obtained in a simpler way based on 
the underlying conformal symmetry of the problem [3] • 

Now, we want to find corrections to ipo- This is nothing 
but a deformation of the minimal surface © due to the 
exponential factor in the metric ([1]). For our purposes, it 
is sufficient to consider ipi and ip2- 

We begin with the leading correction. Using §5§ and 
© , we find that ipi obeys the following linear differential 
equation 



(l-* 2 )<-^i-4i 2 =0. 



(10) 



Again, we look for a solution which is regular at t — 1 
and satisfies "0i(O) = 0. This requirement is met by fixing 
the two free parameters of a general solution (solution to 
(1 - t 2 )^'{ - (2/t)tp' 1 = 0). Finally, takes the form 



(1 - t 2 ) ^'-(^-i^ + ^i)^) -^a-4*Vi = 0, 

(12) 

where ipi is given by (jlip . 

As before, we look for a solution that is regular at t = 1 
and vanishing at t — 0. A simple algebra shows that this 
requirement allows one to fix the two free parameters of 
a general solution. As a result, we have 



ip 2 {t) = 16Li 2 



l+t 



, . 28 
-41n2 t- ~rt 



+ 8 ln(l + t) ln(l + t) + 2 ln(l - t) + t 2 + t - 



- 16 ln21n(l -t) - -tt 2 +8 In 2 2. (13) 
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Here Li 2 is the dilogarithm function. We have used that 
Li 2 (i) = i 7 r 2 -|ln 2 2. 

Having found the minimal surface, we can now cal- 
culate the corresponding area. Like tp, it can also be 
expanded in powers A such that S = J2 n =o ^n^ n - 

To begin with, let us consider the leading contribution 
to S. It follows from © that it is simply given by 



S o (e)=0 



1 dt 
. H 2 



fa 



(14) 



As in 13, there is a divergence in Sq. To make sense 
of the minimal area, we have regularized the integral. In 
doing so, we cut the integral off at the lower bound. In 
terms of r, this means keeping only the part of the surface 
with r < — (f ) 2 - After subtracting the divergent 
term, we find 



S = -a. 



(15) 



Now we would like to go beyond the leading approx- 
imation and consider the corrections. The computation 
is straightforward but somewhat lengthy. One can eval- 
uate Si by reducing the integrals to elementary ones and 
boundary terms. In evaluating the boundary terms, it is 
useful to use the following asymptotics: ipi 
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we get 



t 6 and 

(4 In 2 - 4f )t 6 as t -> 0. At the end of the day, 
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(16) 



ipi(t) =4<-2t 2 -41n(l + f). 



(11) 



Now we move on to the next-to-leading correction. In 
this case equation (fIU|) is replaced by 



Here we set e = because the integrals converge. Note 
that K w 0.14. 

Having derived the leading and sub-leading correc- 
tions, we can now write down the expectation value of 
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the circular Wilson loop of radius a. In the approxima- 
tion we are using, the result is given by 

W ~ exp jg ^1 - ^ca 2 - ^kc 2 ^ + 0{a 6 )j J . (17) 

Qualitative features. Evaluation of the gluon conden- 
sate in terms of small Wilson loops is common on the 
lattice, see @, Q . Namely, one starts with the smallest 
Wilson line possible, that is plaquette and represents it 
as follows: 

\nW = -J2cna , s-^ZG 2 s 2 + 0( S 3 ), (18) 

n 

where the sum means the perturbative contribution and 
the gluon condensate introduced in l| is given by the 
coefficient in front of s 2 (s — ira 2 ). Moreover, we put 
the number of colors N c — 3 and the factor Z is de- 
fined in terms of the perturbative beta function as Z^ 1 = 
0/0W(g), with /3« the one-loop beta function. 

Compare representations (fT8|) and (fT"7f for the expec- 
tation value of the Wilson loop in the dual and direct 
formulations of the Yang-Mills theories. There are terms 
proportional to s 2 in the both cases and this allows for a 
unique determination of the gluon condensate, see below. 

There is a subtle point: a dimension two correction 
that scales like s. It is obviously present in (fT7| but some- 
what disguised in (fT5|) . where it does appear due to the 
UV renormalons [l4j . Phenomenologically the quadratic 
correction is known to be related to the confinement ef- 
fects [13]. Unlike the gluon condensate, however, this 
correction is not related to a matrix element of dimension 
two operator. The dual formulation reveals the univer- 
sality of the quadratic correction, relating it to form of 
the metric in the fifth dimension. This is to be consid- 
ered as a phenomenological success of the gauge/string 
duality. 

The dual formulation settles also another mystery 
of the lattice determination of the gluon condensate. 
Namely, usually the Aq C£) term, or condensate is related 
to an infrared renormalon, or factorial divergence of the 
perturbative series in (|18p . However, the first 16 terms 
of the perturbative expansion evaluated in Ref. Q do 
not exhibit any factorial divergence while allowing for a 
consistent determination of the gluon condensate. In the 
dual formulation, the gluon condensate is associated with 
small values of the fifth coordinate z, see above. Thus, 
the dual formulation avoids the problem of the infrared 
renormalon and provides a qualitative explanation for the 
properties of the perturbative expansion in (fT5)) . 

Determining the gluon condensate. Now, the gluon 
condensate can be read off from (fT7|) . and so 

G 2 = — a k%Z- 1 c 2 . (19) 



Making an estimate requires some numerics. First, 
as in [9], the value of c is fixed from the slope of the 
Regge trajectory of p(n) mesons. This gives c = 0.91 ± 
0.011 GeV 2 if the resonances with n > 3 are used. At this 
point the error is estimated to be less than 2%. Next, the 
overall constant g is adjusted to fit the slope of the lin- 
ear term of the Cornell potential [ll[ . This assumes that 
(eg/47r)c = l/o 2 , where l/o 2 is the slope of the Cornell 
potential. As known, it is a fitting parameter defined 
from the quarkonia spectra. Using the values of [IH , it is 
estimated to be 1/a 2 = 0.186 ± 0.021 GeV 2 . Putting all 
this together, we find g = 0.94 ±0.11. Here, the error is 
now of order 11%. At first glance, 1 — Z _1 can be made 
arbitrary small by letting the scale a go to 0. This is 
however in conflict with the fit of 1/a 2 done within the 
range of 0.2 — 1 fm. What saves the day is the fact that 
the linear term in the quark potential is valid at least 
up to 0.05 fm [16]. Then, the size of the leading cor- 
rection (which is RS independent) syjj is 0.08 — 0.14 
depending on the number of quarks. In our estimate of 
the condensate we will assume that Z~ x ~ 1 within 20% 
accuracy. 

Assembling all the factors, we have 

G 2 = 0.010 ± 0.0023 GeV 4 . (20) 

This is our main result. It is surprisingly close to the orig- 
inal phenomenological estimate 0.012 GeV 4 of [H, though 
somewhat smaller than another phenomenological esti- 
mate 0.024 GeV 4 of [17]. For comparison, the lattice cal- 
culations give bigger values like 0.04 GeV 4 @. 

Conclusions. Using AdS/QCD correspondence, we 
found corrections of dimension two and four to the lead- 
ing zero-point-fluctuations contribution to small Wilson 
loops. The estimated value of the gluon condensate is re- 
markably close to the phenomenological estimates. The 
dual formulation provides also means for universal de- 
scription of the quadratic correction which has never been 
found in terms of the direct formulation. Our results 
suggest that a dual formulation might be very good for 
describing physics at strong coupling. 
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